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Abstract 

In this paper, a ne-w method of Hoo observer design for Lipschitz nonlinear systems is pro- 
posed in the form of an LMI optimization problem. The proposed observer has guaranteed 
^ I decay rate (exponential convergence) and is robust against unknown exogenous disturbance. 
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In addition, thanks to the linearity of the proposed LMIs in the admissible Lipschitz con- 
stant, it can be maximized via LMI optimization. This adds an extra important feature 
to the observer, robustness against nonlinear uncertainty. Explicit bound on the tolerable 



^\ , nonlinear uncertainty is derived. The ne-w LMI formulation also allo-ws optimizations over 

I the disturbance attenuation level {Hao cost). Then, the admissible Lipschitz constant and 

the disturbance attenuation level of the Hoo filter are simultaneously optimized through 



LMI multiobjective optimization. 



Keywords: Lipschitz nonlinear systems, Robust observers, Nonlinear H^o filtering, LMI 
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5^ , 1 Introduction 

The design of nonlinear state observers has been an area of constant research for the last three 
decades and as a result, a -wide variety of design techniques for nonlinear observers exist in the 
literature. Despite important progress, many outstanding problems still remain unsolved. A 
class of nonlinear systems of special attention is the so-called Lipschitz systems in -which the 
mathematical model of the system satisfies a Lipschitz continuity condition. Many practical 
systems satisfy the Lipschitz condition, at least locally. Roughly speaking, in these systems, 
the rate of gro-wth of the trajectories is bounded by the rate of gro-wth of the states. Observer 
design for Lipschitz systems -was first considered by Thau in his seminal paper [1] where he 
obtained a sufficient condition to ensure asymptotic stability of the observer. Thau's condition 
provides a very useful analysis tool but does not address the fundamental design problem. 
Encouraged by Thau's result, several authors studied observer design for Lipschitz systems 
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[21 O HI [5l [6]. All these methods share a common structure for the error dynamics of the 
nonlinear systems; namely the error dynamics can be represented as a linear system with a 
sector bounded nonlinearity in feedback. This type of problems are both theoretically and 
numerically tractable because they can be formulated as convex optimization problems [7], 
[S]. Raghavan formulated a procedure to tackle the design problem. His algorithm is based 
on solving an algebraic Riccati equation to obtain the static observer gain [2]. Unfortunately, 
Raghavan's algorithm often fails to succeed even when the usual observability assumptions 
are satisfied. Raghavan showed that the observer design might still be tractable using state 
transformations. Another shortcoming of his algorithm is that it does not provide insight into 
what conditions must be satisfied by the observer gain to ensure stability. A rather complete 
solution of these problems was later presented by Rajamani [3]. Rajamani obtained necessary 
and sufficient conditions on the observer matrix that ensure asymptotic stability of the observer 
error and formulated a design procedure, based on the use of a gradient based optimization 
method. He also discussed the equivalence between the stability condition and the minimization 
of the Hoo norm of a system in the standard form. However, he pointed out that the design 
problem is not solvable as a standard i^oo optimization problem since the regularity assumptions 
required in the Hoc framework are not satisfied. Using Riccati based approach, Pertew et. al. 
[6] showed that the condition introduced in [3] is related to a modified Hoc norm minimization 
problem satisfying all of the regularity assumptions. It is worth mentioning that the H^o 
problem in [3] is associated with the nominal stability of the observer error dynamics while no 
disturbance attenuation is considered. Moreover, in all of the above references, the system model 
is assumed to be perfectly known with no uncertainty or disturbance. In order to guarantee 
robustness against unknown exogenous disturbance, the nonlinear i^oo filtering was introduced 
by De Souza et. al. [9l |T0] via the Riccati approach. In an Hqo observer, the /32-induced gain 
from the norm-bounded exogenous disturbance signals to the observer error is guaranteed to 
be below a prescribed level. On the other hand, the restrictive regularity assumptions in the 
Riccati approach can be relaxed using linear matrix inequalities (LMIs). In this paper, we 
introduce a novel nonlinear i?oo observer design method for Lipschitz nonlinear systems based 
on the LMI framework. Our solution follows the same approach as the original problem of Thau 
and proposes a natural way to tackle the problem, directly. Unlike the methods of [21 [3l [6], 
the proposed LMIs can be efficiently solved using commercially available software without any 
tuning parameters. In all aforementioned references, the Lipschitz constant of the system is 
assumed to be known and fixed. In this paper, the resulting LMIs are formulated such that 
to be linear in the Lipschitz constant of the nonlinear system. This adds an important extra 
feature to the observer, robustness against nonlinear uncertainty. Maximizing the admissible 
Lipschitz constant, the observer can tolerate some nonlinear uncertainty for which an explicit 
norm-wise bound is derived. In addition to this robustness, we will extend our result such that 
the observer disturbance attenuation level (the Hoo feature of the observer) can be optimized 
as well. Then, both the admissible Lipschitz constant and the disturbance attenuation level are 
optimized simultaneously through multiobjective convex optimization. The rest of the paper is 
organized as follows: Section 2, introduces the problem and some background. In Section 3, the 
LMI formulation of the problem and our observer design algorithm are proposed. The observer 
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guaranteed decay rate and robustness against nonlinear uncertainty are discussed. In Section 

4, we expand the result of Section 3, to an Hoo nonlinear observer design method. Section 

5, is devoted to the simulators optimization of the observer features through multiobjective 
optimization. In section 6, the proposed observer performance is shown in some illustrative 
examples. 

2 Preliminaries and Problem Statement 

Consider the following continuous-time nonlinear system 

x{t) = Ax{t) + u) Ae M"^" (1) 
y{t) = Cx{t) C G (2) 

where x G M"',u G W^,y G W and ^{x,u) contains nonlinearities of second order or higher. 
We assume that the system ([1|-(I2|) is locally Lipschitz in a region T> including the origin with 
respect to x, uniformly in u, i.e.: 

-$(x2,n*)K7||xi -X2II Vxi(A:),X2(A;) gP (3) 

where ||.|| is the induced 2-norm, u* is any admissible control signal and 7 > is called the 
Lipschitz constant. If the nonlinear function <I> satisfies the Lipschitz continuity condition 
globally in M", then the results will be valid globally. Consider now an observer of the following 
form 

S:{t) = Ax{t) + ^{x,u) + L{y -Cx). (4) 

The observer error dynamics is given by 

e{t) ^ x{t) - x{t) (5) 
e{t) = {A- LC)e{t) + $(x, u) - $(x, u). (6) 

The goal is to find a gain, L, such that: 

• In the absence of disturbance, the observer error dynamics is asymptotically stable i.e.: 
limj^oo e(t) = 0. 

• In the presentence of unknown exogenous disturbance, a disturbance attenuation level is 
guaranteed. {Hoo performance). 

The result is simple and yet efficient with no regularity assumption. The observer error 
dynamics is asymptotically stable with guaranteed decay rate (the convergence is actually 
exponential as we will see). In addition, the observer is robust against nonlinear uncertainty and 
exogenous disturbance. The dismissible Lipschitz constant which as will be shown, determines 
the robustness margin against nonlinear uncertainty, and the disturbance attenuation level (the 
Hao cost), are optimized through LMI optimization. 
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3 An Algorithm for Nonlinear Observer Design 



In this section an LMI approach for the nonlinear observer design problem introduced in Section 
2 is proposed and some performance measures of the observer are optimized. 

3.1 Maximizing the Admissible Lipschitz Constant 

We want to maximizes the admissible Lipschitz constant of the nonlinear system (l)-(2) for 
which the observer error dynamics is asymptotically stable. The following theorem states the 
main result of this section. 



Theorem 1. Consider the Lipschitz nonlinear system ([2P-(0j along with the observer 
The observer error dynamics ^ is (globally) asymptotically stable with maximum admissible 
Lipschitz constant if there exist scalers e > and > and matrices P > and F such that 
the following LML optimization problem has a solution. 



s.t. 



min{^) 

A^P + PA- C'^F'^ - FC < -I - el 



p ki 



once the problem is solved 



> 



L = P^^F 



(7) 
(8) 



(9) 

7' = max(7) = (10) 
Proof: Suppose Q = F The original problem as discussed in section 2, can be written as 



s.t. 



min(Amaz(^')) 

{A - LCfP + P^{A - LC) = -I 

l-2-i.\^ax{P) > 

P > 



(11) 

(12) 
(13) 



which is a nonlinear optimization problem, hard to solve if not impossible. We proceed by 
converting it into an LMI form. A sufficient condition for existence of a solution for (llip is 

3 e > 0, {A- LCfP + P^{A - LC) < -I - el. (14) 
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The above can be written as 

A^P + PA- C^L^P - PLC <-I-eI (15) 
which is a bihnear matrix inequaUty. Defining the new variable 

F = PL^ L'^P^ = L^P = (16) 

it becomes 

A^P + PA- C'^F'^ - FC < -I - el (17) 
In addition, since P is positive definite (j{P) = Xmax{P)- So, from (fT2]) we have 

-a{P) < i- (18) 

which is equivalent to 

{—fl-P^P>0 (19) 
27 

using Schur's complement lemma 

I J-T P 1 

> (20) 



fl P 
P 



defining ^ = ^, ([8]) is achieved. A 

Proposition!. Suppose the actual Lipschitz constant of the system is 7 and the maximum 
admissible Lipschitz constant achieved by Theorem 2, is 7* . Then, the observer designed based 
on Theorem 2, can tolerate any additive Lipschitz nonlinear uncertainty with Lipschitz constant 
less than or equal 7* — 7. 

Proof: Assume a nonlinear uncertainty as follows 

^A(a;,u) = <I>(x,u) + A<I>(x,u) (21) 
x{t) = Ax{t) + ^/:^{x,u) (22) 

where 

||A$(xi,u) - A$(x2,n)|| ^ A7||2;i - X2||. (23) 

Based on Schwartz inequality, we have 

\\^a{.xi,u) - ^/s.{x2,u)\\ < \\^{xi,u) - $(x2,n)|[ + ||A$(xi,n) - A$(x2,'u)|| 

< 7||xi - X2II + A7IIX1 - X2II. (24) 
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According to the Theorem 1, ^/\{x,u) can be any Lipschitz nonhnear function with Lipschitz 
constant less than or equal to 7*, 

||^>A(a;i,n) - <^a{x2,u)\\ < 7*lki - X2II (25) 

so, there must be 

7 + A7 < 7* ^ A7 < 7* - 7. A (26) 

Remark 1. If one wants to design an observer for a given system with known Lipschitz 
constant, then the LMI optimization problem can be reduced to an LMI feasibility problem 
(just satisfying the constraints) which is easier. 

^From Theorem 2, it is clear that the gain L obtained via solving the LMI optimization 
problem, can lead to stable error dynamics for every member in the class of the Lipschitz non- 
linear functions with Lipschitz constant less than or equal to 7*. Thus, it neglects the structure 
of the given nonlinear function. It is possible to take advantage of the structure of the <I>(x,ti) 
in addition to the fact that its Lipschitz constant is 7. According to Proposition 1, the margin 
of robustness against nonlinear uncertainty is 7* — 7. The Lipschitz constant of the systems can 
be reduced using appropriate coordinates transformations. The transformation matrices that 
are picked are problem specific and they reflect the structure of the given nonlinearity [2]. The 
robustness margin can then be modified through coordinates transformations. Finding the Lip- 
schitz constant of a function is itself a global optimization problem, since the Lipschitz constant 
is the supremum of the magnitudes of directional derivatives of the function as shown in [TT] 
and [12] . If the analytical form of the nonlinear function and its derivatives are known explicitly, 
any appropriate global optimization method may be applied to find the Lipschitz constant. If 
only the function values can be evaluated, a stochastic random search and probability density 
function fitting method may be used [13| . 

3.2 Guaranteed Decay Rate 

The decay rate of the system ([6]) is defined to be the largest /? > such that 

limexp(/3t)||e(t)|| =0 (27) 

holds for all trajectories e. We can use the quadratic Lyapunov function V{e) = e^Pe to es- 
tablish a lower bound on the decay rate of the (0). If ^ ~2/3y(e(t)) for all trajectories, 
then V{e{t)) < exp(-2/3i)y(e(0)), so that ||e(t)|| ^ exp(-/?t)«;(P)5 ||e(0)|| for ah trajectories, 
where k{P) is the condition number of P and therefore the decay rate of the ([6]) is at least /3, 
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[8]. In fact, decay rate is a measure of observer speed of convergence. 



Theorem 3. Consider Lipschitz nonlinear system along with the observer 

The observer error dynamics (0j is (globally) asymptotically stable with maximum admissible 
Lipschitz constant and guaranteed decay rate (3, if there exist a fixed scaler /? > 0, scalers e > 
and > and matrices P > and F such that the following LMI optimization problem has a 
solution. 

min{£) 

s.t. 

A^P + PA + 2/3P - C'^F^ - FC < -I - el (28) 



p ki 



> (29) 



once the problem is solved 



L = P^^F (30) 
7* = max(7) = (31) 

Proof: Consider the following Lyapunov function candidate 

V{t) = e^{t)Pe{t) (32) 

then 

V{t) = f{t)Pe{t) + e^{t)Pe{t) = -e^Qe + 2e^P($(x, u) - $(x, u))^ . (33) 
To have V{t) ^ — 2/3y(t) it suffices ([33]) to be less than zero, where: 

{A - LCfP + P'^iA - LC) + 2I3P = -Q. (34) 

The rest of the proof is the same as the proof of Theorem 2. A 

4 Robust Hoo Nonlinear Observer 

In this section we extend the result of the previous section into a new nonlinear robust Hoo 
observer design method. Consider the system 

x{t) = Ax{t) + ^{x,u) + Bw{t) (35) 
yit) = Cx{t) + Dw{t) (36) 



7 



where 'w{t) G i22[0, oo) is an unknown exogenous disturbance, suppose that 

z{t) = He{t) 



(37) 



stands for the controhed output for error state where H \s a, known matrix. Our purpose is to 
design the observer parameter L such that the observer error dynamics is asymptoticahy stable 
and the following specified i?oo norm upper bound is simultaneously guaranteed. 

||-z|| < /"ll^^ll- (38) 

The following theorem introduces a new method for nonlinear robust i?oo observer design, we 
first present an inequality that will be used in the proof of our result. 

Lemma 1 [14^ . For any x,y ^ R" and any positive definite matrix P G M"^", we have 

2x^y < x'^Px + y^P'^y (39) 

Theorem 4- Consider the Lipschitz nonlinear system (^^-^J^ with given Lipschitz con- 
stant 7, along with the observer The observer error dynamics is (globally) asymptotically 
stable with decay rate (3 and minimum 2^2{w e) gain, /j,, if there exist fixed scaler (3 > 0, 
scalers a > 1, e > and C > and matrices P > and F such that the following LMI opti- 
mization problem has a solution. 



min{C,) 



s.t 



A' P + PA + 2/3P -C'F' -FC < -al - el 

> 



P 

H'^H + + - - PB-FD 



Once the problem is solved 



L = P-^F 
IJ* - min(/i) = 



< 



Proof: The observer error dynamics will be 

e{t) = {A- LC)e{t) + ^{x, u) - $(x, u) + {B - LD)w 



(40) 
(41) 

(42) 



(43) 
(44) 

(45) 
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consider the following Lyapunov function candidate 

V{t) = e^{t)Pe{t) (46) 

then 

V{t) = e^{t)Pe{t) + e^{t)Pe{t) = -e^Qe 
+ 2e^P{^{x, u) - $(x, u)f + e'^iPB - FD)w + nF {B'^ P - D'^ F^)e (47) 



where, Q is as in (|34p . We select Q = al. li w = the error dynamics is as Theorem 2, so the 
LMIs d?]) and ([8]) which for Q = al will become 

A^P + PA + 2(3P - C'^F^ - FC < -al - el (48) 

-/ P 



> (49) 



27 
P 

are sufficient for the asymptotic stability of the error dynamics. Having a > 1, (jlSp always 
implies (09]). 

Based on Rayleigh inequality 

e^Qe < Xmax{Q)e^e (50) 

Using Lemma 1 we can write 

2e^P{<^>{x, u) - <^ix, u)) < e^Pe + ($(x, u) - ^>(x, u)f PP-^P{<^>{x, u) - «>(x, u)) 

= e^Pe + {^{x, u) - ^{x, u))'^P{^{x, u) - ^{x, u)) (51) 
based on Rayleigh inequality we have 

||e^Pe|l < XmaxiP)\\ef = A„a,(P)e^e (52) 

||($(x,n) -^>(x,n))^P($(x,n) -$(x,n))|| < X^acciP)\\'^ix,u) - <l>ix,u)f 

< l^Xmax{P)\\ef = -f^Xn,ax{P)e^e (53) 
therefore, from the above and (1181). 



2e^Pi<^{x,u)-^x,u)) < {l + j^)X^^,{p)e^e<li^ + -)e^e. (54) 

According to ()50p and (j54p and knowing that Q = al, we have 

V{t) < -(7 + - - 2a)e^e + e^{PB - FD)w + w^{B'^P - D^F^)e. (55) 
2 7 



Now, we define 

poo 

J= {z^ z - C.'uFw)dt 
Jo 

therefore 

poo 

J < {z'^z - Cw^w + V)dt 
Jo 

it follows that a sufficient condition for J < is that 



VtG[0,cx)), z' z-Cw^w + V <0 



but we have 



z'^z - (w'^w + V < e^H^He - (w^w + Ve^H^He + -(7 + ^ - 2a)e^e 

2 7 



+ e' {PB - FD)w + uF{B'^P - D^F^)e - (w' w 



(56) 
(57) 

(58) 



e 
w 



H + + ^ - 2a)I PB-FD 



BTp - D^F^ 



-CI 



e 
w 



(59) 



Thus, a sufficient condition for J < is that the above matrix which is the same as ()42p be 
negative definite. Then 

z^z - (w'^w <0 ^ \\z\\ < y^\\w\\ (60) 

Up until now, we have the LMIs ([48]), ([20]) and (|42]). If these LMIs are all feasible, then 
the problem is solvable and the observer synthesis is complete. However, (j20p can be slightly 
modified to improve its feasibility. We proceed as follows: 

Inequality ()5ip can be rewritten as follows 

2e^P($(x, u) - $(x, u)) < 2jXmax{P)e^e (61) 
following the same steps, the matrix in ()59p will become 

H'^H+[2-fXmax{P)-a]I PB-FD 



< 0. 



(62) 



The above matrix can not be used together with (j48p and (j49p because it includes P as one of 
the LMI variables, thus resulting in a problem that is not linear in P. It can, however, give us 
another insight about Xmax{P)- According to the Schur's complement lemma, ([62]) is equivalent 
to 

-CI <0 (63) 
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H^H + [2-f\rnax{P) - a]I + ^{PB- FD){PB - FDf < 0. (64) 
The third term in the above is always nonnegative, so it is necessary to have 

H''H+ [2-fX,nax {P)-a]I<0 (65) 
but as for any other symmetric matrix, for H^H, we have 

\^in{H^H)I < H^H < X^ax{H^H)I (66) 
or according to the definition of singular values 

a^iH)I < H'^H < a^{H)I (67) 
therefore, a sufficient condition for ([65]) is 

a\H) + 2^X^axiP)-a<0 (68) 

or 

\max{P) < (69) 

but (fTSl) must be also satisfied. To have both (fTHI) and (i69l). it is sufficient that 



27 



which is equivalent to (j4ip . A 



Remark 2. Similar to Remark 1, if one wants to design an observer for a given system with 
known Lipschitz constant and with a prespecified fi, the LMI optimization problem is reduced 
to an LMI feasibility problem. 

Remark 3. As an additional opportunity, we can first maximize the admissible Lipschitz 
constant using Theorem 3, and then minimize ^ for the maximized 7, using Theorem 4. In this 
case, according to Proposition 1, robustness against nonlinear uncertainty is also guaranteed. 
In the next section, we will show that how 7 and fi can be simultaneously optimized using 
convex multiobjective optimization. It is clear that if no decay rate is specified, then the term 
2/3P will be eliminated from LMI ()40p in Theorem 4. 
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5 Combined Performance using Multiobjective Optimization 



The LMIs proposed in Theorem 4 are hnear m both admissible Lipschitz constant and dis- 
turbance attenuation level and as mentioned earlier, each can be optimized. A more realistic 
problem is to choose the observer gain matrix by combining these two performance measures. 
This leads to a Pareto multiobjective optimization in which the optimal point is a trade-ofF 
between two or more linearly combined optimality criterions. Having a fixed decay rate, the 
optimization is over 7 (maximization) and fi (minimization), simultaneously. The following 
theorem is in fact a generalization of the results of [21 |3l IU O El |15] , and [9] (for systems of 
class dl])-©) in which the Lipschitz constant is assumed to be known and fixed and the result 
of [7] in which a special class of sector nonlinearities is considered. 



Theorem 5. Consider the Lipschitz nonlinear system. (3^- (3^ along with the observer 
The observer error dynamics is (globally) asymptotically stable with decay rate f3 and si- 
multaneously maximized admissible Lipschitz constant 7* and minimized £,2{w — t- e) gain, /i*, 
if there exist fixed scalers < A < 1 and (3 > 0, scalers a>l, e>0, ^>0 and C > and 
matrices P > and F such that the following LMI optimization problem has a solution. 



min [A • C + (1 - A)C] 



s.t 



A^P + PA + 2/3P - C^F^ -FC < -al - el 



p 



> 



H'^ H + - 2a)I L PB-FD 



-2CL 






-a 



< 



Once the problem is solved, 



L = P-^F 



7 



* A 

/i = mm 



max(7) = ^ ^ 

(^) = x/C 



(71) 
(72) 

(73) 



(74) 
(75) 
(76) 



Proof: The above is a scalarization of a multiobjective optimization with two optimality 
criteria. Since each of these optimization problems is convex, the scalarized problem is also 
convex [16] . The rest of the proof is the same as the proof of Theorem 4 where the LMI (j73p 
is obtained from the LMI (j42p using the Schur's complement lemma. A 
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6 Illustrative Examples 



In this section the high performance of the proposed observer is shown via three design examples. 



Example 1. Consider the fohowing observable (A,C) pair 



A 



1 

1 -1 



1 



The result of the iterative algorithm proposed in ^ is 



L 



0.49 

69.5523 11.5679 



while using our proposed method in Theorem 2, 



L 



1.1933 
56.8334 21.9074 



which means that the admissible Lipschitz constant is improved by a factor of 2.42. 



Example 2. The following system is the unforced forth-order model of a flexible joint 
robotic arm as presented in [2], [5], [3]. The reason we have chosen this example is that it is an 
important industrial application and has been widely used as a benchmark system to evaluate 
the performance of the observers designed for Lipschitz nonlinear systems. 



X 





-48.6 


19.5 



1 

-1.25 







48.6 

1 

-19.5 



x + 







-3.33 sin(x3) 



10 
10 



X. 



The system is globally Lipschitz with Lipschitz constant 7 = 3.33. Noticing the structure 

of $ that has a zero entry in three of its channels, Raghavan [2], proposed the coordinates 
transformation x = Tx, where 

T = diag [1,1, 4:,0.1] (77) 

under which the transformed system has Lipschitz constant 7 = 0.083. Using Theorem 3, 
7* = 0.4472 in the original coordinates and 7* = 2.4177 in the transformed coordinates. The 
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observer gain L, is obtained in the transformed coordinates and computed in the original 
coordinates as L = T^^L. Assuming 



B 

D 
H 



0.2 

1111 

0.1 0.25 

0.5/4x4, 



T 



using Theorem 4 we get, /i* = 0.5753, a = 2.0517, e = 0.0609, and finally the observer gain will 
be 



L 



33.4865 129.9249 59.89713 108.2134 
38.5694 282.8603 102.1561 171.0910 



Figure [H shows the true and estimated values of states. The actual states are shown along 
with the estimates obtained using Raghavan's and Aboky's methods and our proposed LMI 

1 T 



optimization method. The initial conditions for the system are x(0) 

1 T 







-10 2 



and 



those of the all observers are x(0) 



As seen in Figure [H the observer 
las the best convergence of the three. 



1 -0.5 

designed using the proposed LMI optimization method 
Note that in addition to the better convergence, the proposed observer is an H^o filter with 
maximized disturbance attenuation level while the observers designed based on the methods of 
[21 El m El E] can only guarantee stability of the observer error. 



Example 3. In this example we show the usage of of the multiobjective optimization of 
Theorem 5 in the design of //qo observers. Consider the following system 



X 



X 



y 



Xi X2 

1 
-1 -1 



X 



-6xf - 6xf X2 - 2x| - 2xf 



1 



X. 



(75) 



The systems is locally Lipschitz. Its Lipschitz constant is region-based. Suppose we consider 
the region V as follows 



V 



(xi,X2) G I xi < 



0.25| 
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2 
1.5 
1 

0.5 


-0.5 





Actual 

LMI 

Raghavan 

- — Aboky 


ft \^ 
}' 





X 



2 3 

time(Sec) 




2 3 

time(Sec) 



Figure 1: The true and estimated states of Example 2 

in which the Lipschitz constant is 7 = 0.4167. We choose 

H = 0.5/ 

r 1 T 
B= 11 

D = 0.2 
/3 = 0.05 
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and solve the multiobjective optimization problem of Theorem 5 with A = 0.9. We get 



7* = 0.5525 
/i* = 1.1705 
a = 1.6260 



L 



2.2435 X 10" 



23.7025 13.7272 



The true and estimated values of states are shown in Figure [2j We have assumed that 

1 T 



x(0) 
wit) 



-0.2 -1.45 

1 ^ 

0.25 -2 
0.15exp(-t)sin(t). 




Figure 2: The true and estimated states of Example 3 in the presence of disturbance 

The values of 7*, fi* , norm of the observer gain matrix, (t{L), and the optimal trade-off 
curve between 7* and /i* over the range of A when the decay rate is fixed (/3 = 0.05) are shown 
in Figure [3l The optimal surfaces of 7* and /i* over the range of A when the decay rate is 
variable are shown in Figures S] and [5l respectively. 
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Figure 3: 7*, fi* and o"(-L), and the optimal trade-off curve with /3 = 0.05 

7 Conclusions 

A new method of robust observer design for Lipschitz nonlinear systems proposed based on 
LMI optimization. The Lipschitz constant of the nonlinear system can be maximized so that 
the observer error dynamics not only be asymptotically stable but also the observer can tolerate 
some additive nonlinear uncertainty. In addition, the result extended to a robust Hoo nonlin- 
ear observer. The obtained observer has three features, simultaneously. Asymptotic stability, 
robustness against nonlinear uncertainty and minimized guaranteed Hoo cost. Thanks to the 
linearity of the proposed LMIs in both admissible Lipschitz constant and the disturbance atten- 
uation level, they can be simultaneously optimized through convex multiobjective optimization. 
The observer high performance showed through design examples. 
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